Local symmetry and magnetic anisotropy in multiferroic MnWC>4 and 
antiferromagnetic C0WO4 studied by soft x-ray absorption spectroscopy 
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Soft x-ray absorption experiments on the transition metal Z/2,3 edge of multiferroic M11WO4 and 
antiferromagnetic C0WO4 are presented. The observed linear polarization dependence, analyzed by 
full-multiplet calculations, is used to determine the ground state wave function of the magnetic Mn 2+ 
and Co 2+ ions. The impact of the local structure and the spin-orbit coupling on the orbital moment 
is discussed in terms of the single-ion anisotropy. It is shown that the orbital moment in C0WO4 is 
responsible for the collinear antiferromagnetism, while the small size of spin-orbit coupling effects 
make spiral magnetic order in MnWC>4 possible, enabling the material to be multiferroic. 

PACS numbers: 71.20.Be, 78.70.Dm, 71.70.Ch, 75.30.Gw 



The combination of magnetic order and ferroelectric- 
ity in multiferroic compounds is a topic that has gained 
high interest in modern research.— Magnetic transitions 
directly linked to a spontaneous electrical polarization 
have been investigated in detail for materials like i£Mn03 
(R=Y, Tb, Gd, Dy) 2 .^ and TbMn 2 0^. Among the mul- 
tiferroic manganates, M11WO4 takes a special place as it 
does not contain a rare-earth ion, leaving Mn formally as 
the only magnetic ion£i£ 

Multifcrrocity in these materials is explained by the 
spin current model where the electrical polarization is 
caused by a spiral order of the magnetic moments. The 
spin rotation axis e does not coincide with the magnetic 
propagation vector Q, leading to a non- vanishing spon- 
taneous polarization P oc e x Q. Here we investigate the 
conditions for the occurrence of spiral magnetic order re- 
sponsible for the rich physics of M11WO4. We performed 
soft x-ray absorption spectroscopy (XAS) on M11WO4 to 
determine the local electronic and magnetic properties of 
the Mn ions. The results are also compared to measure- 
ments on the isostructural compound C0WO4, which is 
a collinear antifcrromagnet and not a multiferroic mate- 
rial. We analyze the role of the spin-orbit coupling and 
its impact on the magnetic structure. The magnitude 
of the single-ion anisotropy is calculated using the full- 
multiplet configuration-interaction approach. This forms 
the key to understand why in M14WO4 spiral magnetic 
order is possible and why it is suppressed in C0WO4. 



I. EXPERIMENTAL 

Single crystals of M11WO4 and C0WO4 were grown 
from the melt. During the single crystal growth process 
of M11WO4 special care was taken to keep manganese 
in divalent state. This was achieved by the use of high 



growth temperatures and the renouncement of melt sol- 
vents. Using the top seeded growth technique and low 
cooling rate (0.08 K/h), ruby- red transparent single crys- 
tals were obtained starting from a growth temperature of 
1574 Kj 1 ^ Single crystals of C0WO4 were obtained from 
Na2W207 melt solution with a small surplus of WO3 
(molar ratio C0WO4 : Na 2 W 2 07 : W0 3 = 1:2: 0.5), 
starting at 1363 K and applying a cooling rate of 3 K/h. 

The XAS spectra were collected at the Dragon Beam- 
line of the National Synchrotron Radiation Research 
Center in Hsinchu, Taiwan. The samples were cleaved 
in-situ in an ultra-high vacuum chamber with pressures 
in the 10~ 10 mbar range, guaranteeing the high sur- 
face quality required to perform bulk-rcprcscntative XAS 
studies in the total electron yield mode. The degree of 
linear polarization was 99 ± 1%, with an energy resolu- 
tion of approximately 0.3 eV. The spectra were collected 
at room temperature in the paramagnetic phase of both 
compounds. Single crystals of MnO and CoO were mea- 
sured simultaneously in a separate chamber as references 
for M11WO4 and C0WO4, respectively. 



II. SPECTROSCOPIC RESULTS 

Fig. [T] shows the room temperature Mn-L 2i 3 XAS spec- 
tra of M11WO4 taken with the E vector of the light paral- 
lel to the a, b, and c crystallographic axes. The spectrum 
of a MnO single crystal is also included for reference pur- 
poses. The spectra are dominated by the Mn 2p core-hole 
spin-orbit coupling which splits the spectrum roughly in 
two parts, namely the L3 {hv ps 638-645 eV) and L 2 
(hv as 649-654 eV) white line regions. The line shape 
strongly depends on the multiplet structure given by the 
Mn 3d-3d and 2p-Zd Coulomb and exchange interactions, 
as well as by the local crystal fields and the hybridization 
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FIG. 1: (color online) Top panel: experimental and theoreti- 
cal Mn Z/2,3 XAS spectra of MnWCUwith the E vector of the 
light parallel to the a, b, and c crystallographic axes. The 
spectrum of MnO is included as reference. Bottom panels: a 
close-up revealing the polarization dependence of the spectra. 



with the O 2p ligands. Unique to soft XAS is that the 
dipole selection rules are very sensitive in determining 
which of the 2p 5 3d n+1 final states can be reached and 
with what intensity, starting from a particular 2p 6 3d n 
initial state (n = 5 for Mn 2+ ) u i 12 . This makes the tech- 
nique extremely sensitive to the symmetry of the initial 
state, i.e., the valence and the crystal field states of the 
ions. 

Comparing the MnW04 spectra with the MnO, one 
can immediately observe that the spectral features are 
similar with the L3 main peaks at identical energies. This 
indicates directly that the Mn ions in MnW04 are also 
in the high-spin Mn 2+ electronic configuration as the Mn 
ions in MnO. Yet, the clear differences between the spec- 
tra of MnW04 and MnO, as seen e.g. at the low and 
high energy shoulders of the L3 white line, provide a 
hint to a crystal field level scheme different from Oh for 
the MnW04 system. 

A small but clear polarization dependence for MnW04 
was found, as can be seen in the bottom panels of Fig. [TJ 
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FIG. 2: (color online) Experimental and theoretical Co £2,3 
XAS spectra of C0WO4. The spectrum of multi-domain CoO 
single crystal is included as reference. 

where a close-up is shown of the spectra taken with the 
three different polarizations. The intensities as well as 
the energy positions of several peaks vary with the po- 
larization of the incoming light. This polarization depen- 
dence is yet another indication that the local symmetry 
of the Mn ions in MnW04 is lower than Oh ■ 

Fig. [2] depicts the Co-L 2j3 XAS spectra of C0WO4 
taken with the E vector of the light parallel to the a, 
b, and c crystallographic axes, together with the spec- 
trum of a multi-domain CoO single crystal as reference. 
Following the same argumentation as given above for 
the Mn £2,3 edge, from the strong similarities between 
the C0WO4 and the CoO spectra we can conclude that 
C0WO4 contains Co 2+ ions with a 3c? 7 high-spin config- 
uration like in CoO. In comparing the C0WO4 with the 
M11WO4, one observes that the polarization dependence 
is much larger. In the following we will discuss the dif- 
ferent origins of the linear polarization dependence for 
the two compounds and its implications for their very 
different magnetic properties. 



III. LOCAL ELECTRONIC STRUCTURE 

To interpret and understand the spectral line shapes 
and their polarization dependence, we have performed 
simulations of the atomic-like 2p 6 3d n 2p 5 3d n+1 (n = 5 
for Mn 2+ and n = 7 for Co 2+ ) transitions using the 
well-proven configuration-interaction cluster modelAi~— 
Within this method we have treated the Mn or Co ion 
within an MnOe or CoOg cluster, respectively, which in- 
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FIG. 3: (color online) The structure of MnW0 4 and C0WO4. 
For the calculations of the XAS spectra, a local coordinate 
system (x, y, z) was introduced, in which the axes point 
roughly along the Mn/Co-0 bonds. The global crystallo- 
graphic coordinates system can be transformed into the lo- 
cal coordinate system by rotating along b by ~ 34° and then 
along the c axis by 45° (see text for more details). 



eludes the full atomic multiplet theory and the local ef- 
fects of the solid. It accounts for the intra-atomic 3d-3d 
and 2p-3d Coulomb interactions, the atomic 2p and 3d 
spin-orbit couplings, the local crystal field, and the O 2p- 
Mn 3c? or O 2p-Co 3d hybridization. This hybridization 
is taken into account by adding the 3d' 1+1 L and 3cP +2 L 2 
etc. states to the starting 3d™ configuration, where L 
denotes a hole in the O p ligands. 

To facilitate the setup and the interpretation of the 
calculations, a local coordinate system was introduced, 
with its axes {x, y, z) pointing roughly along the Mn-0 
and Co-0 bonds. To transform the global coordinates 
(a, b, c) into the local system, first a rotation by 33.9° 
along b is applied, resulting in a new system with (a', 
c'). This system is then rotated by 45° along the old 
axis c, transforming (a', 6', c') into (x, y, z). The local 
coordinate system is shown in Fig. [31 The coordinates 
of the oxygen ions as well as the Mn-0 and Co-0 dis- 
tances are summarized in Table HI The crystal structure 
has been taken from Ref. [H]. There are two distorted 
octahedra in the unit cell and they are connected to each 
other by inversion, thereby contributing to the polariza- 
tion dependence in the same manner. 

The Mn-0 and Co-0 bond lengths were used to 
estimate the hybridization strength using Harrison's 
description^ resulting in values for pda and pdir as 
shown in Table Q] The Slater-Kostcr formalism- 6 * 
provides the angular dependence of the hybridization 
strength. Values for the crystal fields were tuned to 
find the best match to the experimental spectra. Pa- 
rameters for the multipole part of the Coulomb inter- 
actions were given by the Hartree-Fock values^ while 
the monopolc parts (JJdd, U p d) were estimated from pho- 
tocmission experiments on MnQii and previous work on 



MnWO-4 






X 


y 


z 


Mn-0 dist. 


pda 


pdir 


1 


0.25 


0.24 


2.15 


2.18 


-1.18 


0.54 


2 


0.25 


0.24 


-2.15 


2.18 


-1.18 


0.54 
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0.31 


-2.07 


-0.02 


2.09 


-1.36 


0.63 
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-2.06 


0.36 


0.02 


2.09 


-1.36 


0.63 
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0.31 


2.24 


-0.08 


2.27 


-1.02 


0.47 
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2.25 


0.25 


0.08 


2.27 


-1.02 


0.47 



C0W04 






X 


y 


z 


Co-0 dist. 


pda 


pdn 


1 


0.19 


0.19 


2.10 


2.12 


-1.07 


0.49 
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0.19 


0.19 


-2.10 


2.12 


-1.07 


0.49 
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0.24 


-2.02 


-0.04 


2.03 


-1.24 


0.57 
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-2.02 


0.25 


0.04 


2.03 


-1.24 


0.57 


5 


0.24 


2.14 


-0.03 


2.16 


-1.01 


0.47 


6 


2.14 


0.23 


0.03 


2.16 


-1.01 


0.47 



TABLE I: Coordination of the Mn 2+ and Co 2+ ions in their 
distorted oxygen octahedra. All oxygen positions are given 
in A in the local coordinate system specified in fig. |3j where 
the numbering of the atoms is also defined. The hybridization 
strengths pda and pdir are in units of eV. 



CoO.— The simulations were carried out using the pro- 
gram XTLS 8.3,— and the parameters used are listed in 
Refs.tilM 

The calculated polarization-dependent spectra for the 
Mn L 2i 3 edge of MnW04 are plotted in Fig. [TJ One can 
observe that the general line shape of the experimental 
spectra is very well reproduced. Equally important, also 
the small energy shifts in the peak positions as well as the 
small variations in the peak intensities as a function of 
polarization can all be simulated. This indicates that we 
have been able to capture the local electronic structure 
of the Mn ion in MnW04 with great accuracy. We now 
will look into this in more detail. 

The ground state of a Mn 2+ ion in Oh symmetry is 
6 ^4i, and this will not split or change upon lowering the 
crystal symmetry. The electronic charge distribution of 
this high-spin half-filled shell ion remains spherical. The 
presence of lower symmetry crystal fields will then show 
up most clearly when an electron is added (or removed) 
from the valence shell. To illustrate such a situation, 
we depict in the left panel of Fig. [4] the total energy 
level diagram of a Mn + (3d 6 ) ion using the crystal and 
ligand field parameters which reproduce best the Mn £2,3 
spectra of MnW04 as described above. One can readily 
observe a splitting of about 0.7 eV between states with 
the added electron in the t2 g vs. the e g orbitals. Upon 
switching off the 3c? spin-orbit interaction, sec Fig. @] (a) , 
one can also identify an intra-i2g splitting of about 65 
and 95 meV as well as an intra-e ff splitting of 50 meV. 
All of the orbitals are of mixed type as a result of the 
low local symmetry. The two orbitals lowest in energy 
are approximately 1/V2(xz + yz) and 1/V2(xz — yz), 
resembling the xz and yz orbitals, rotated by ±45° along 
z. This is a result of the local S2 symmetry axis of the 
MnOg cluster, which lies in the xy plane with an angle 
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FIG. 4: Energy level diagrams for (a,b) Mn+ (3d 6 ) and (c,d) 
Co 2+ clusters, excluding and including the 3d spin-orbit cou- 
pling f. Only the states up to 1 eV are shown. 



of roughly 45° to the oxygen bonds. 

Also in the Mn £2,3 XAS process one extra electron 
is added into the valence shell. The lowest lying peak 
at 639.6 eV, see Fig. [IJ involves the excitation of the 
2p core electron into the t2 g orbitals. Different polariza- 
tions access different orbitals with different probabilities. 
For the 640.6 eV peak, the polarization dependence is 
reflected in terms of shifts in its energy position. If one 
could have measured the polarization dependence along 
the local coordinates of one Mn ion, one may expect to 
observe peak shifts of 65 and 95 meV associated with the 
intra t2 g splittings shown in Fig. |4] Yet, having two Mn 
ions in the unit cell, and measuring along the global a, 6, 
and c crystallographic directions, the shifts between the 
effectively composite peaks become reduced. In the ex- 
periments and in the simulations we are left with shifts 
of about 20 and 30 meV as can be seen in the middle 
bottom panel of Fig. [TJ For the 639.6 eV peak, the shifts 
due to variation of the polarization are more difficult to 
quantify since there is also a variation in the intensity of 
the peaks. The latter is apparently caused by the rather 
complicated multiplct effects involving the 2p core hole, 
as demonstrated by the excellence of the match by the 
simulations, see the left bottom panel of Fig.[TJ These in- 
tensity variations are yet quite small, and integrated over 
the entire £2,3 range, would have been even identical to 
zero if the non-cubic hybridization effects were absent. 

The simulations for the polarization dependent spectra 
for the Co £2,3 edge of C0WO4 are shown in Fig. [2] Like 
in the Mn case, we have been able to achieve a satisfying 
fit to the experimental spectra: all features are well re- 
produced, including the strong polarization dependence. 



We will now analyze the local electronic structure of the 
Co ion on the basis of the parameters used in these sim- 
ulations. 

For the Co 2+ ions, the ground state of the 3c? 7 config- 
uration in Oh symmetry is 4 7i (~ ^2g e g)- Unlike for the 
6 Ai of the Mn, this 4 Ti state will be split upon going to 
lower crystal symmetry. The degeneracy in the orbital 
part, i.e. in the tig sub-shell, will be lifted. This can 
be seen from the total energy level diagram shown in the 
right panel of Fig. 01 The calculation with the 3d spin- 
orbit interaction switched off, see Fig. H](c), reveals the 
presence of three low- lying quartets within the first 0.1 
eV. In an one-electron language, each of them would cor- 
respond to the hole occupying one of the three ^orbitals. 
The intra-£2g splittings are 40 meV and 80 meV. These 
values are quite close to those of the MnWO.4 case re- 
flecting the similar crystal structure. The importance of 
the lifted degeneracy is that the orbital occupation will 
no longer be isotropic, i.e. the Co ion charge distribu- 
tion will be highly non-spherical. This then explains the 
strong polarization dependence in the absorption spectra. 

It is also important to notice how the low-lying states 
of Co 2+ are influenced by the 3d spin-orbit coupling. The 
energy level diagram in Fig. |4] (d) reveals that there is 
indeed a large amount of mixing due to spin-orbit inter- 
action. The spin-orbit coupling constant of C, = 66 meV 
is of the same order of magnitude as the splittings in 
the tig subshcll. All this has direct consequences on the 
magnetism, as will be discussed in the next section. 



IV. SINGLE-ION ANISOTROPY 

From the results of the simulation, we find that the 3c? 
spin-orbit interaction has a major influence on the en- 
ergies and nature of the low lying states of the C0WO4 
system. The magnetic moment will then not only have 
a spin contribution but also an appreciable orbital one. 
This results in a strong coupling of the direction of mag- 
netic moment to the crystal structure. We can make an 
estimate of the single ion anisotropy of the Co 2+ ion. 
For that we use the same parameters with which we have 
reproduced the XAS spectra and the polarization depen- 
dence therein, and we apply a small exchange field of 
f-BH ex = 10~ 6 eV (H ex ?s 200 Oe) to align the magnetic 
moment. Here we would like to note that the results 
hardly changes with increasing the strength of the field, 
i.e. they are robust even beyond [IbH £X = 10~ 2 eV. 

Fig. [5] shows the expectation value of the magnetiza- 
tion parallel to the applied exchange field My as a three- 
dimensional representation and in a polar plot, the latter 
in local coordinates. It can be seen that the single-ion 
anisotropy is of easy-axis type, with the axis lying along 
4> = 135° at 8 = 90° in the local coordinate system. In 
the direction of the easy axis, the moment is completely 
directed along the magnetic field, and the perpendicu- 
lar component of the magnetization vanishes. The easy 
axis is also perpendicular to the S2 symmetry axis of the 
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FIG. 5: (color online) Calculated magnetic anisotropy of Co 2+ 
in C0WO4. (a) Three-dimensional representation of the mag- 
netic anisotropy in the unit cell. The surface depicts the ex- 
pectation value of the magnetic moment My pointing paral- 
lel to a small exchange field applied to the ion. (b) Two- 
dimensional cut of My in the plane of easy and intermediate 
magnetic axis (corresponding approximately to the xy plane 
in local coordinates), (c) Cut through the plane of hard and 
intermediate magnetic axis. 



distorted octahedron, reflecting the orbital character of 
the two lower lying ligand field levels. The total magnetic 
moments of the hard and easy axes are 1.2 (1b and 3.0 hb, 
respectively, with an orbital contribution of 0.1 fiB and 
0.9 [is- The easy axis of the total moment coincides 
directly with that of the orbital moment, which is also 
pointed out earlier by Brunoi^I It has also been reported 
that for low symmetries, the orbital and total moment 
could differ in their behavior due to the magnetic dipolc 
moment T z *£ In the case of C0WO4, however, with the 
symmetry of the CoOg groups being not too far from cu- 
bic, this T z is around —0.01 eV, and does not influence 
the anisotropy significantly. 

In C0WO4, the single- ion anisotropy is strong enough 
to lead to an effective pinning of the spins in the mag- 
netically ordered phase. Isotropic superexchange will be 
the dominating magnetic interaction, although in addi- 
tion anisotropic exchange is possible, as the bonds lack 
inversion centers. But the Dzyaloshinskii-Moriya (DM) 



term D- (Si x S2), that favors canted magnetic moments, 
is small compared to the energy involved in the magnetic 
anisotropy of C0WO4, reflected by the large difference 
in moment of more than 1 ^b- The superexchange will 
then minimize the energy in the magnetically ordered 
phase by aligning the moments along the easy axis. The 
easy axis is identical for both CoOg clusters in the unit 
cell, and thus the single-ion anisotropy itself will not lead 
to canting between the magnetic moments of the Co 2+ 
ions. The result is a collincar antifcrromagnet with an 
easy axis approximately in the ac plane ss 40° off from a. 
This direction is illustrated in Fig. [5] and is close to the 
one found in the experiment by neutron diffraction^ 

Returning to MJ1WO4 case: to the leading order, the 
local spherical 6 Ai ground state does not produce an or- 
bital moment. The spin is then free to point in any di- 
rection. This is reproduced in the calculation, where the 
response of the magnetic moment to an exchange field is 
practically isotropic. Yet, in second order, the spin-orbit 
coupling does affect the Mn-0 hopping. This produces 
a non-vanishing moment perpendicular to the exchange 
field, albeit in the order of 10 -2 /ib, being much lower 
than in the case of the spin-orbit active C0WO4. Con- 
sequently, the DM term is on the same energy scale as 
the anisotropy, and spin canting can easily occur. The 
propagation direction of the spiral magnetic order is in 
the local xy plane, meaning that the second order pro- 
cess leads to a hard axis along z. This demonstrates the 
importance of the single-ion anisotropy in this system. 
The form of the anisotropy, however, depends strongly 
on the magnetic exchange. This explains that even three 
different magnetically ordered phases exist in the mate- 
rial at low temperatures^ Additionally to the hard axis, 
the easy plane also has two principal axes for the mag- 
netic anisotropy with a small difference between the two. 
This, in combination with the subtle dependency of the 
anisotropy on the exchange field, makes incommensura- 
bilities likely to appear. 



V. CONCLUSION 

The single-ion anisotropics extracted from our soft x- 
ray absorption data and the corresponding full-multiplet 
calculations provide a natural explanation for the mag- 
netic ordering phenomena occurring in C0WO4 and 
MnW04. In C0WO4, the spin-orbit coupling and a 
non- vanishing orbital moment overwhelm the anisotropic 
Dzyaloshinskii-Moriya exchange and cause the formation 
of a collincar ordering of the moments. Contrary to that, 
although being isostructural, MnW04 shows only a small 
magnetic anisotropy which is comparable in energy to 
the Dzyaloshinskii-Moriya interactions. This will cant 
the spins in such a way that the spiral propagates per- 
pendicular to the direction of smallest magnetic moment, 
allowing a modulated spin spiral which causes the ferro- 
electricity in the material. 
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